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ABSTRACT 



We calculate the GeV spectra of GRB afterglows produced by inverse Comp- 
ton scattering of the sub-MeV emission of these objects. We improve on earlier 
treatments by using refined afterglow parameters and new model developments 
motivated by recent Swift observations. We present time-dependent GeV spec- 
tra for standard, constant parameter models, as well as for models with energy 
injection and with time-varying parameters, for a range of burst parameters. We 
evaluate the limiting redshift to which such afterglows can be detected by the 
GLAST LAT, as well as AGILE. 

Subject headings: gamma rays: bursts - radiation mechanisms: non-thermal 



Introduction 



GRB afterglow obser vations from X- rays to radio are generally well described by an 
external shock model 



(e.g. 



Meszaros 2006, for a recent review). However, model fits to the 



data still leave uncertainties in some of the parameters of the basic external shock model, 
and the bolometric luminosity depends on the p oorly known GeV range spectrum. The 
Gamma Ray Large Area Space Telescope GLAST ( IMcEnery fc GLAST Mission Teamll2006l ) 
is expected to be launched at the end of 2007. The Large Area Telescope (LAT) covers the 
energy range from 20 MeV to 300 GeV, while the Gamma Ray Burst Monitor (GBM) covers 
the range from 8 keV to 25 MeV. The effective area of the LAT is about 7 times larger than 
of the previous EGRET experiment at GeV energies. AGILl£ was successfully launched on 
April 23, 2007, with an energy range of 30 MeV to 50 GeV (ITavani et al.ll2006l ). It is hoped 
that the higher photon statistics at GeV energies of GLAST and AGILE will lead tighter 
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constraints on the burst parameters, and an improved understanding of the GeV afterglow 
spectra. 

GLAST may also be able to test recent developments in the understanding of GRB after- 



glows, motivated by observations with th e Swift satellite (e.g. lNousek et al.ll2006l ; IZhang et al. 



20061 ; iFox fc Meszarosll2006l ; IZhangj|2007l ). One such development is the presence of a shallow 
decay phase o f the X-ray afterg low, which may be due to refreshed shocks or late energy 
injection (e.g. IZhang et al.l 120061). or al ternatively, it may be due to a change in time of the 



Ioka et al.ll2006l ). We investigate here the effects of such new features 



shock parameters (e.g. 
on the expected GeV spectrum. 

Another question of great interest is how far can GLAST detect the MeV to GeV emis- 
sion from such bursts, both in the basic model and the case where new features such as 
the above are present. This requires a detailed calculation of the GeV spectrum as a func- 
tion of time, with allowance for the changes in the dynamics implied by the injection, time 
variability, etc. The most conservative and widely cons idered mechanism for producing pho 



tons i n this range is inverse Compton (IC) scattering (IPanaitescu fc Meszaroslll 998; 



btani 



19981; 


Wei & Lu 19981: Chiansr & Dermer! 19991 


Panaitescu & Kumar! 2000: ' 


Sari & Esin 2001 




Zhansr & Meszarosl2001: Wane et al.koOl. 20061: Wei & Fan 20071: Fan et al 


2007: Galli & Piro 


2007' 


. Another potential mechanism is hadronic cascades following proton acceleration 


(Bottcher & Dermerlll998; Zhang & Meszaros 


2001: Frarile et al.ll2004: Gupta & Zhanell2007 





This mechanism is less well constrained, depending on the efficiency of proton acceleration; 
it may be important in the prompt phase, but its parameter regime is small and generally 



outside of the afterglow parameter fit range (IZhang fc Meszarosll200ll ). so it is not considered 



he re. The maximum d i stanc e to which GLAST could detected GRB afterglow was discussed 



by lZhang fc Meszarod (120011 ) for the basic standard model, using a simplified analytical bro- 
ken power-law approximation to the IC spectrum. This resulted in an IC-to-synchrotron 
peak-flux ratio which is overestimated by a factor ~ 10, compared to a more accurate cal- 
culation, as we discuss in this paper. The usefulness of this ratio is that it allows simple 
predictions for the detectability in the GeV range based on the lower energy measurements. 
Here we discuss the validity of the simple analytical approximations, compared to more 
accurate numerical calculations of afterglow spectra at GeV energies. 

In § [2] we describe the afterglow synchrotron-inverse Compton model used for the nu- 
merical computations, a comparison between numerical and analytical approximations being 
given in the Appendix. In § [3] we then present both numerical IC spectra and their appro- 
priate analytical approximations, for the basic GRB afterglow model and for the extended 
models including new Swift-motivated elements such as energy injection or time-varying 
parameters, and evaluate their detectability as a function of redshift for the GLAST LAT 
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instrument and for AGILE. 



2. Afterglow Synchrotron-inverse Compton Spectra at GeV Energies 



The afterglow of a GRB, due to the external shock as it slows down in the external 
medium, produces synchrotron radiation in the X-ray to MeV range, w 



Comp ton upscattered into the GeV-TeV range (IMeszaros et al.l Il994 



199 8|). More specific calc u lation of the IC G e V ran ge are, e.g. those of 



rich is then 


inverse- 


Bottcher & Dermer 


Sari & Esin 


( 


2001 





Zhang fc Meszarod (120011 ). IPe'er fc Waxmanl (120041 ). etc. We describe the afterglow models, 
as usual, by the total isotropic energy E 52 ,iso = -E/10 52 erg (for the case of energy injection see 
below), the external density n, a jet opening half-angle 9, electron equipartition parameter 
e e , magnetic equipartition parameter e# and electron energy index p. The other parameter 
of relevance in synchrotron-IC models is the scattering Y-parameter which is defined as the 
luminosity ratio of IC to synchrotron, usually given by 



y = (-l + v/4e e /e B + l)/2 



The initial Lorentz factor Tq of the burst is not needed as a parameter, since in the asymptotic 
blast wave regime the Lorentz factor follows from the scaling law, 



I\t) = (17E/102ATrnm p c b t 



,5j.3\1/8 



(2) 



Our numerical calc ulations of the sp ectra and fluence curves use the basic synchrotron-IC 
equations given in iGou et al.l (120071 ). extending now to the GeV range. We also consider 
in this range the spectral effects of the photon-photon opacity effe cts, which impose cutoffs 
depending on the spe ctrum and density of lower energy photons (IBaring fc Harding! 119971 ; 
Lithwick fc SarilhoOlh . 



We calculate the synchrotron-IC spectra of three different types of GRB afterglow mod- 
els, and evaluate their detectability with the GLAST Large Area Telescope (LAT) and wit h 
AGILE. These calculations improve on previous calculations, e.g. (IZhang fc Meszarosll20 01). 
in sey eral respects. First, in the "standard" afterglow model such as used in (jZhang &: Meszaros 
20011 ). we use the exact IC spectrum instead of the broken power-law approximation, and 
the peak flux ratio is taken as Fl instead of F2 (see Appendix). Second, we include the 
complete spectral regimes, not only the commonly used v a < (u m , u c ) regimes where u a , u m , 
and v r are the s ynchrotron self-absorption, injection, and cooling frequencies, respectively 
( IGou et al.ll2007l ). This assures that the GRB afterglows evolve through the correct regimes 
at all times. Third, we consider Swift-motivated additions to the standard afterglow model, 
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such as a continued energy injection model, and a varying afterglow parameter model (mo- 
tivated by the presenc e of a shallow decay phase, and a high apparent radiative efficiency, 
Mfezarosll2006 for a review). 



see e.g. 



The details of the three models calculated are as follows. (A) A standard afterglow 
model, with all par ameters remaini ng constant during its evolution (for a detailed description 
of this model, see Gou et al .120071 ) . (B) A continuous energy injection model, which is a 
widely considered model to explain the shallow decay phase commonly seen in Swift light 
curves. For this, we assume that the total kinetic energy increases over time with a power- 
law index, E oc t 1 " 9 , before the break time t = 10 4 seconds and the break time here is 
defined as the one when the shal low decay phase ends. Fits to Swift observations indicate a 
value q ~ 0.5 (jZhang et al.ll2006l ). (C) An evolving parameter model, which is an alternative 
model for explaining the shallow decay phase, based on the assumption that the electron 
equipartition parameter e e increases with the time as oc t a ( lloka et al.ll2006l ) before the break 
time, as for the energy injection model. For all three models, we assume that they have the 
same parameters at late times, i.e. after the break time. Since the flux has to be integrated 
over the observation time, we set the observation time to be one half of the final time of 
observation since the trigger (e.g. if the observational data stop at t — 10 5 seconds, the 
integration time is from i = 5x 10 4 seconds to t — 10 5 seconds). This is consistent with the 
GLAST observation characteristics, as well as those of AGILE, in the point-source observing 
mode, where due to earth occultation, only about 50% of the orbit time is used for the burst 
observation. 

To determine the limiting redsh ift to which a burs t can b e detected, we calculate the 
instrumental fluence threshold as in IZhang fc Meszarod (120011 ). using the instrument char- 
acteristics of the GLAST LAT and AGILE. For a flux sensitivity $ m ph cm -2 s _1 over an 
exposure time T and a point source observed over an effective observation time t e ff (in 
seconds) in an energy band centered around a photon energy E, the fluence threshold is esti- 
mated as F t hr ~ [&o(T /t e ff Y^ 2 ]Et e ff where [Q (T /t e ff) 1 / 2 ] is the sensitivity for the effective 
observation time t e // because the sensitivity scales as y/t e ff for the longer observations where 
the sensitivity is limited by the background. Due to the occultation by the earth, the effective 
observation time t e ff is normally < 50% of the total orbit time, t b. s , for both GLAST and 
AGILE (or equivalent to the observation time after the burst), namely, t e ff = ^t b s where 
the observing efficiency rj is taken to be 77 = 0.5. Hereafter, unless otherwise stated, the small 
"£" without the subscript refers to the observation time t b s for simplicity. For GLAST we 
use the fluence threshold listed in the updated instrument performance documents^]. Consid- 
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ering the integral sensitivity above 100 MeV for GLAST LAT to be ~ 4 x 10 _9 ph cm -2 s _1 
for an effective observation time of one year in the sky-survey mode, the fluence threshold 
is 1.0 x 10-H 1 / 2 erg cm for the long-time observation in the sky-survey mode. For GRB 
afterglows, in most cases GLAST will perform a pointed observation rather than the survey 
mode observation. In this mode, GLAST will keep the GRB position always at the center of 
the LAT field of view for as long as possible, and this will improve the sensitivity by a factor 
of 3-5 (depending on where the GRB lies with respect to the orbital plane; an object which 
lies at the orbit pole will not be occulted by the Earth and can be continuously observed; 
J. McEnery 2007, private communication). Therefore, taking the improvement factor of 3, 
the fluence threshold for a GRB observation is 3.4 x 10 _9 t 1//2 erg cm -2 for the long-time 
observation. The short-time fluence threshold can be defined by the criterion that at least 5 
photons are collected and it depends on the effective area of the instrument which is around 
6000 cm 2 at 400 MeV for GLAST LAT, so it is 5.3 x 10~ 7 erg cm" 2 (the transition time when 
the short-time sensitivit y and long-time sensitivit y meet is 2.4 x 10 4 seconds). Compared to 



the previous estimate of IZhang fc Meszaroa (120011 ) for the GLAST sensitivity, the short-time 



observation sensitivity here is roughly the same, but the long-time sensitivity has changed 
from 1.2 x lO" 9 ^ 2 to 3.4 x lO" 9 ^ 2 erg cm , increased by a factor of 3. 

The energy range of AGILE is somewhat lower than that of GLAST LAT. Its flux 
sensitivity above 100 MeV is ~ 3 x 10~ 7 ph cm" 2 s" 1 for a point-source observation over 



an effective period of 10 6 seconds ( ITavani et al.l 120061 ). Thus, taking the observing effi- 



ciency i] = 0.5, the fluence threshold can be estimated as 3.0 x 10" 7 (T/teffY^Eteff 
1.3 x 10" 7 t 1/ ' 2 ergs cm -2 at an average energy of 400 MeV, where T = 10 6 seconds is the 
exposure time corresponding to the given sensitivity. The fluence threshold is for the long- 
time observation. The shorter-time fluence threshold can be obtained similarly as above, 
~ (5/550) x 400 MeV ~ 5.8 x 10~ 6 ergs cm -2 where we have taken the effective area for 
AGILE to be 550 cm 2 . In summary, for AGILE, the fluence threshold for the long-time 
observation (i.e., t > 1870 seconds) is 1.3 x 10~ 7 t 1//2 ergs cm -2 and 5.8 x 10~ 6 ergs cm -2 for 
the short-time observation, the transition time being ~ 1870 seconds. 



3. Detectability of GRB Afterglows with GLAST and AGILE 



The initial nominal set of parameters for the stan dard model (A) used here are the same 
as for the standard model of IZhang fc Meszarod (120011 ): p = 2.2, e e = 0.5, = 0.01, E 52 ,iso = 
l,n = l cm -3 . An additional feature is that we also assume a jet opening half-angle 9 = 0.14, 
which does not affect the flux at early times. The other parameters are as for model (A). 
For the injection model (B), the kinetic energy is assumed to increase following the relation 
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E oc t l - q = t - 5 where q is taken to be 0.5 based on the Swift observation fits (jZhang et al.l 



2006). For the model (C) with evolving parameters, we assumed e e to follow the relation e e oc 



t - 5 Jlokaet aDl2006h . other parameters being the same as for model (A). In the alternative 



models (B) and (C), either the kinetic energy or the electron equipartition parameter starts 
out with a smaller value as for (A), but at late times end up with the same values as the 
standard model (A). The transition time at which the energy injection or the e e evolution 
stops is set at t — 10 4 seconds. 

Fig [H shows the results for the three models above, using the nominal set of parameters. 
Panel (a) shows the partial fluence, defined here as the energy flux integrated over the time 
intervals [t — At, t], as a function of t, where the t = £ & s is the observation time counted after 
the trigger, adopting a nominal integration time At = 0.5t throughout. The partial fluence 
curves shown correspond to the three different GRB a fterglow models, (A) the standard 



model, which is the same as the type II GRB model in IZhang fc Meszarod (120011 ); (B) the 
energy injection model; (C) the evolving parameter model. As can be seen in panel (a), for 
a burst at a low redshift z = 0.32 the GeV emission from all three models can be detected 
by GLAST up to a time t ~ 1.5 x 10 5 seconds (the thick solid line indicates the GLAST 
LAT sensitivity). Note that the GeV emission from the standard model is higher than that 
from the two other models. This is because all the models end up with the same energy 
and same parameters at late times, which means the injection starts with lower energy and 
the evolving parameter begins with the lower e e at the beginning. Panel (b) shows the 
synchrotron and IC spectra of the standard model (A) at times t = 10 2 , 10 3 , 10 4 , 10 5 , and 
10 6 seconds. Hereafter, unless otherwise stated, we always calculate the spectra at these time 
epoches. The fluxes around TeV (10 12 eV) show the effects of inclusion of the photon-photon 
absorption within the sources. The upper curves are the flux without the 7-7 absorption, 
and the lower curves are the flux after internal absorpti on. For this we have used t he optical 



depth to internal 7-7 interactions of Equation (20) in IZhang fc Meszarod (120011 ). For the 
relatively low compactness parameters of the external afterglow shock discussed here, the 77 
cutoff becomes important above ~ TeV energies, which is more of interest for ground-based 
air Cherenkov telescope observations than for space detectors. The lower panels (c) and (d) 
show the redshift dependence of the GeV emission for all three models, at t ~ 1.1 x 10 3 
and t ~ 2 x 10 4 seconds. At t ~ 1.1 x 10 3 seconds, the limiting redshift is z ~ 0.4 for the 
standard model and z ~ 0.22 for the other two models. At t ~ 2 x 10 , the limiting redshift 
is around z ~ 0.45 for all the models. 

Note that while the usual fluence is defined as flux integrated over the observation time 
since the trigger, which always increase with time, the partial fluences shown in panel (a) 
first increase and eventually decrease. This is because the afterglow flux decreases with time 
t, and for the partial fluences the integration time starts at 0.5t and ends at t. This is done 
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to check the limiting redshift to which afterglows can be detected for typical observations 
at different epochs t with some uniform criterion for the integration time. The snapshot at 
the time 2 x 10 4 seconds lies where the partial fluence is roughly flat in time, during which 
period the limiting redshift reaches its maximum (although the partial fluence within the 
flat phase varies by a factor < 2, the limiting redshift changes only slightly). Other snapshot 
epochs were chosen around one decade earlier or later than the typical maximum redshift 
epoch. 

In Fig [2] we consider an alternative set of parameters. The motivation for this is that 
the parameters of the standard model shown in Figure [fl E^ 2 ,iso = 1, and e e = 0.5, differ 



somewhat from the 'statistical average' y a 



ues q uoted for low- redshift GRBs, E 52) j e t ~ 0.1 



and e e ~ 0.1, e.g. iPanaitescu &: Kumar! (120011 ) . In order to check the sensitivity of the 



detectability of GRBs to variations in these parameters, we perform the same calculation 
suing the values p = 2.2, e e = 0.2, = 0.01, E§ 2> i so = 10,n = 1 cm -3 , and 9 = 0.14, the 
results being shown in Fig ([2]). For these 'average' parameters, the limiting redshift is z ~ 0.8 
for all three GRB models at t ~ 2.0 x 10 4 seconds. 

In Fig [3] we show the corresponding results for AGILE. This is a smaller-scale mission 
than GLAST, launched in April 23, 2007, and it is interesting to compare its detectability 
limits with those of GLAST. AGILE has a different energy range (30 MeV to 50 GeV) and 
has a relatively narrower energy band than the GLAST LAT (20 MeV to 300 GeV), as well 
as a lower effective area. Thus the observed fluxes and partial ffuences are expected to be 
lower for AGILE in contrast to GLAST. This is seen in Panel (a) of Fig [3J The dashed 
line is for AGILE, and the solid line is for GLAST, showing that it is hard for AGILE to 
detect a burst with the typical parameters at z = 0.32, while GLAST could detect it until 
around 2 days. Panels (b), (c) and (d) of Fig [3] show the detectability with AGILE and with 
GLAST at different time epochs. At t = 1.1 x 10 3 seconds, AGILE can detect bursts up 
to z ~ 0.25, and GLAST can detect bursts to z ~ 0.8. At the time t = 2.0 x 10 4 seconds, 
the limiting redshifts are 0.15 for AGILE, and 0.8 for GLAST(same as Fig [2]), respectively. 
At t — 1.4 x 10 5 seconds, the limiting redshift for AGILE is apparently well below z = 0.1, 
while the limiting redshift for GLAST can still reach up to 0.5. We see that the limiting 
redshift for AGILE drops relatively quickly with increasing observation times; the short-time 
sensitivity for AGILE lasts around 10 3 seconds, and after that the sensitivity drops quickly. 
For GLAST, the short-time sensitivity lasts a longer time, ~ 10 4 seconds, and since the 
afterglow GeV flux doesn't change much for up to one day after the trigger, we don't expect 
the sensitivity drop to have much of an effect on the limiting redshift for GLAST. 

In Figure H] we probe the sensitivity of the detectabilty on the total kinetic energy of 
the burst, taking as an example the results for a value of E^ 2 ^ so = 100. This is in the range 
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Fig. I. — Panel (a) The partial fluence curves (defined as flux integrated between 0.5t and t) 
as a function of the observation time t since the trigger, for the three GRB afterglow models 
in the GLAST LAT energy band, at a redshift z = 0.32. Thin solid curves: standard model 
(A) (constant energy); dashed curves: "energy injection" model (B); dot-dashed curves: 
"evolving parameter" model (C). The sensitivity of the GLAST LAT is shown by the thick 
broken solid curve. The parameters for the standard model are E 52 ,iso = 1, e e = 0.5, 
e B = 0.01, p = 2.2,9 = 0.14 and the break time when the shallow decay phase ends is at 
t = 10 4 seconds. Before the break time, e e oc t 5 in the evolving parameter model, and 
the kinetic energy E oc t 05 in the energy injection model. After the break time, those 
parameters will be the same as the ones in the standard model. Panel (b) The synchrotron 
and IC spectrum for the standard model at different time epochs: 10 2 , 10 3 , 10 4 , 10 5 , and 10 6 
seconds, respectively. Above photon energies ~ 10 12 eV, the upper spectral curve represents 
the flux without 7-7 absorption, and the lower curve is with inclusion of this absorption. 
Panel (c) The redshift dependence of the partial fluence for the three models above, evaluated 
at the observation time t ~ 1.1 x 10 3 seconds. The thick solid line is the GLAST sensitivity 
for an integration time 550 seconds (the observing efficiency 0.5 has been taken). The 
intersection of the partial fluence and the sensitivity curve gives the limiting redshift, to 
which the bursts can be detected by GLAST for this integration time, which is is z ~ 0.4 
for the standard model, and z ~ 0.22 for the other two models. Panel (d) The redshift 
dependence of the partial fluence for the three models above, evaluated at t ~ 2.0 x 10 4 
seconds. The thick solid line is the GLAST sensitivity for an integration time 10 4 seconds. 
The intersection gives a limiting redshift for detection of z ~ 0.45 for all three models. 
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Fig. 2. — Panel (a) The partial fluence curves of the same three GRB afterglow models 
over the GLAST energy band, for a different set of parameters generally assumed to be 
'typical', namely kinetic energy E^ 2 ,iso — 10, e e = 0.2 and other parameters the same as in 
Fig ([1]). Panel (b) The redshift dependence of the partial fluence at t = 2 x 10 4 s for the 
same models, giving a limiting redshift z ~ 0.8. Panel (c) The synchrotron-IC spectra for 
the energy injection case. Panel (d) The synchrotron-IC spectra for the evolving parameter 
case. 
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Redshift Redshift 

Fig. 3. — Comparison of GLAST and AGILE detectability, for the standard model with 
the parameters of Fig. [2j Panel (a) The partial ffuence curves for the standard model over 
the energy bands of AGILE (thin dashed line) and GLAST (thin solid line), compared to 
the respective instrument sensitivities (thick dashed line for AGILE and thick solid line for 
GLAST). The parameters are the same as in Fig El Panel (b) The redshift dependence of 
the partial ffuence at t — l.le4 seconds. The limiting redshifts are z = 0.25 for AGILE, 
and z = 0.8 for GLAST, respectively. Panel (c) The redshift dependence of the partial 
ffuence at t — 2.0 x 10 4 seconds. The limiting redshift points are ^ = 0.15 for AGILE, and 
z = 0.72 for GLAST, respectively. Panel (d) The redshift dependence of the partial ffuence 
at if: = 1.4 x 10 5 seconds. The limiting redshift for AGILE is below z ~ 0.1, and the limiting 
redshift for GLAST is around z = 0.5, for this model. 
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of values derived for objects such as GRB 990123 and GRB 050904, which may be called 
hyper-energetic GRB. For the "standard" model case (A) with this energy (see upper panels 
of FigH|) we see that the limiting redshift for a GLAST detection has increased from z ~ 0.8 
to a value z ~ 2.0, assuming that the other parameters remain the same as in Figure EJ 
Thus, hyper-energetic bursts such as GRB 990123, at the observed redshift z = 1.6, should 
be detected by GLAST in the GeV band, if they have the above conventional parameters. 
The other hyper-energetic object, GRB 050904, has a similar kinetic energy as GRB 990123, 
but it was at the much higher redshift z = 6.29, which appears out of range for GLAST. 
GRB 050904 had the most complete set of observational data so far, covering from the BAT 
band, through X-ray, to the optical/NIR and to the radio band. Thus a lot of effort has been 



inves ted in obtaining the best fitting parameters for this burst f Frail et a 



20071 ). Taking the best fitting parameters from model (B) of iGou et al.l (120071 ). our results 



2006; Gou et al. 



here indicate that such GRB 050904-like bursts could be detected by GLAST up to z < 1.0. 
There are two reasons for this relatively modest limiting redshift detectability by GLAST in 
this case: (1) The electron equipartition parameter derived is small, e e = 0.026, which means 
only a fraction of the kinetic energy is radiated; (2) The Compton parameter parameter in 
the fast cooling case is relatively small, Y ~ 1.7, which means that the energy lost via IC 
scattering is comparable to the energy lost via synchrotron radiation. 

In Figure [5] we illustrate the sensitivity of the GeV detectability on the value of the 
Compton Y parameter, showing the fluxes for two values, Y = 2.7 (thin solid line) and 
Y = 6.6 (thin dashed dotted line) in the fast cooling case. The parameters are same as the 
ones of the standard case (A) in Fig. [21 except for the electron equipartition parameter e e , 
which is e e = 0.1 for the Y = 2.7 case and e e = 0.5 for the Y = 6.6 case. One expects a 
higher flux over the GLAST energy band for the Y = 6.6 case because a larger Compton 
parameter means that more energy goes into the GeV band via the IC process, which can 
be seen from the spectrum on the lower panels. The limiting redshifts for the observation 
time t = 2.0 x 10 4 seconds are z = 0.55 and z = 1.2, respectively, for Y = 2.7 and Y = 6.6. 



4. Discussion 

We have calculated the time-dependent GeV synchrotron and inverse Compton spectra 
for three generic types of GRB afterglows which, at lower photon energies, have been used 
to interpret observations from the Swift satellite and other ground-based facilities. These 
models include a standard, constant parameter afterglow model (A), a model with late energy 
injection (B) and a model with varying parameters (C). The spectra and the GeV partial 
fluence curves in the GeV range were used to estimate the detectability with GLAST and 
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Fig. 4. — Standard models of type (A), with different values of the parameters. Panel (a) 
The redshift dependence of the GeV partial ffuence for hyper-energetic GRB 990123-like 
objects, at t — 2.0 x 10 4 seconds. The intersection point gives the limiting redshift z ~ 2 
below which GLAST can detect such afterglows in the GeV band. The model parameters 
are E^,iso = 100, e e = 0.2 and the other parameters are the same as in Fig ([I]). Panel (b): 
The redshift dependence of the pa rtial fluence for model (B) of GRB 050904 using the best 
fitting for this bursts parameters (IGou et al.l 120071 ). p = 2.194, e e = 0.026, eg = 0.0058, n = 
109,6* = 0.111, E 52>iso = 184.8, at the time t = 2.0 x 10 4 seconds. The limiting redshift for 
such GRB 050904-like objects is z 1. Panel (c) The synchrotron and IC spectrum for 
hyper-energetic GRB 990123-like objects. Panel (d) The synchrotron and IC spectrum for 
GRB 050904-like GRBs. 
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Fig. 5. — Detectability limits for GLAST for standard models such as (A) but with different 
Compton Y parameters, Y = 2.7 and Y = 6.6. Panel (a) The partial ffuence curves for 
Y = 6.6 (thin dashed dotted line) and Y = 2.7 (thin solid line), compared to the GLAST 
sensitivity (thick broken solid line). Panel (b) The limiting redshift for an observation time 
t = 2 x 10 4 s and the two Y parameters, giving limiting redshifts z = 0.55, and z = 1.2, 
respectively. Panel (c) The synchrotron and IC spectra for the Compton parameter Y = 2.7 
case. Panel (d) The synchrotron and IC spectra for the Compton parameter Y = 6.6 case. 
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AGILE of bursts in these model categories, for various sets of parameters, at various epochs 
after the trigger and for various observation durations. These model spectra improve on 
previous calculations in several respects. In particular, in the past mainly constant parameter 
spectral models such as (A) were computed; here we have extended these calculations over a 
broader range of input parameters, based on more recent information and statistics. Models 
(B) and (C) have not been previously investigated in the GeV range, and are motivated 
by recent Swift results. T he GeV spe c tra di scussed here were computed numerically, using 



the formalism described in iGou et al.l ( 20071 ). These are compared with previous analytical 



synchrotron-IC spectra of type (A) in the appendix. 

The detectability depends most obviously on the total burst energy E 52 ,iso and on the 
observation time t and integration time At, in addition to the other parameters such as e e , 
€b etc. E.g., for bursts of the standard constant parameters (A) with the nominal values of 
the parameters E^,iso = 1, e e = 0.5, e# = 0.01, the limiting detection redshift with GLAST 
for times t = 1.1 X 10 3 s and 2 x 10 4 s and integration time At = 0. 5t are roughly the same, 
z ~ 0.4 (Fig. [TJ panels b and c). For models with energy injection (B) or varying parameters 
(C), where the final values of Ei SO or e e reach the same value as for model (A) at a later 
time, the detection threshold is somewhat lower for the shorter observation time, as seen in 
the same figure, panel (c), since they start out weaker and build up to be comparable to 
(A) at later times. However, for the longer of the two observation times above, panel d, the 
limiting redshifts are the same. For the more standard values E 52 ,iso = 10, e e = 0.2 (Fig. [2} 
panel b) the limiting redshift for the constant parameter model (A) at t = 2 x 10 4 s goes up 
to z ~ 0.8. 

For AGILE the limiting redshifts are lower than for GLAST due to its lower effective 
area. For a standard burst model (A) with average parameter values E 52tiso — 10, e e = 0.2, 
the limiting redshift is z ~ 0.25 at an earlier observation time t — 1.1 x 10 3 s (Fig. [3J panel 
b), and z ~ 0.15 at t = 2 x 10 4 s. With AGILE bursts can only be detected at relatively 
early times, since the short-time sensitivity for AGILE (where the sensitivity curve is flat 
in time) only lasts, e.g. around ~ 10 3 seconds, versus ~ 10 4 s for GLAST for a burst at 
z ~ 0.3 (Fig. [3l panel a). 

A discrimination between models (A), (B) and (C) based on GeV measurements of the 
spectral evolution in time is possible in principle, as seen e.g. by comparing Fig. 1 panel (c), 
and Fig. 2 panels (c) and (d). However, it will require good energy and time coverage, and 
extensive simulations over a wide range of parameter space, since changes in E iso , Y (e.g. 
Figs. H] and [5]) and the other afterglow parameters needs to be carefully disentangled. 

The limiting redshift naturally increases for larger values of Ei SO , E.g., for GLAST at 
an observation time t = 2 x 10 4 s and a standard model (A) with E§ 2 ^ so = 10 2 , e e = 0.2 it is 
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z ~ 2, while for a for a burst with the parameters of GRB 050904 (high Ei SO but low e e ) it is 
z ~ 1. The limiting redshift also increases with the Compton Y parameter, as illustrated in 
Fig. O which reflects the fact that Y provides a measure of how much energy gets scattered 
into the GeV range. 

Besides the photon-photon absorption inside the source considered in our calculation, 
another interesting absorption process is the photon-photon absorption by the external cos- 
mic infrared background radiation (external absorption), which produces electron/positron 
pairs, and the resulting pairs can IC scatter th e cosmic microw ave background (CMB) pho- 
tons, yielding a delayed MeV-GeV emission (IDai &: Lul 120021 ). It can be shown that the 



external absorption doesn't affect the limiting redshift much. iSalamon fc Steckerl (119981 ) 
show that, at z — 1, the absorption optical depth is r ~ 1 for photons with an energy of 
50 GeV, and r ~ 10 for photons of 300 GeV. If we define the cutoff energy as the photon 
energy corresponding to an optical depth r = 1, the cutoff energy shoul d shift to higher en- 
ergies as one considers lower redshift, based on the numerical results in ISalamon fc Stecker 
(119981 ) . Assuming that (1) the photons above 50 GeV is totally absorbed and there is no 
absorption below 50 GeV; (2) the cutoff energy is roughly the same within the redshift range 
considered; and (3) the flux contribution from the delayed emission is ignored, we can con- 
clude that in this case the flux observed by GLAST will be comparable to that observed by 
AGILE, because the effective observing band will have become similar for both instruments 
(due to external absorption effectively cutting off the GLAST higher energy contribution). 
The limiting redshift for AGILE, as shown in Fig [31 will be approximately the same for the 
case with external absorption as it is without, because absorption is important only above 
the AGILE band. 

The detectability estimates discussed here illustrate the sensitivity to different types 
of model assumptions in current synchrotron-inverse Compton models, when one takes into 
account newer information gleaned from Swift. Calculations using simplified generic mod- 
els sh ow that around tens of Swift-detected GRB per year w ill fall in the LAT field of 
view ( lOmodei &: The GLAST/LAT GRB Science Group! 120061 ) during their prompt emis- 
sion phase. Considering t he Swift-detected, burst redshift distributi on to imply a fraction of 
around 20% below z — 1 (jjakobsson et al. 2006 ; Le Sz Dermer 2006 ). we may roughly expect 
~ 5 Swift-burst prompt detections per year by GLAST. However, since the GeV afterglows 
can last up to a day (e.g. Figs. 1 and 2, panel a), GLAST may actually be able to observe 
more than this number of bursts in the afterglow, as opposed to the prompt phase. De- 
tections with GLAST should test many of the assumptions that go into these models, and 
will provide important new information on the energetics, dynamics and parameters of GRB 
afterglows. 
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5. Appendix 

We discuss here two analytical approximations to the synchrotron-IC spectrum, and 
compare them to the numerically calculated values. The two key elements in the sim- 
ple analytical approximations to synchrotron-IC spectra in the literature are (i) an IC-to- 
synchrotron peak flux ratio, 

p nIC I rsyn /gX 
J max/ J max V / 

expressed, e.g. in erg cm -2 s _1 Hz -1 and evaluated at the frequencies where the synchrotron 
and the IC flux attain their peak value, and (ii) a "Compton parameter" Y, usually taken 
to be 

Y — (—1 + y/4£ e /e B + l)/2 . (4) 

In the GRB literature, the flux ratio ([3]) of the analytical approximations has appeared under 
several forms, two of which are the most relevant for us here. One of these is 



pi = fIC ir 

J max I J i 



esyn 
' max 



( 14/45) a T Rn 



T 



(5) 



(see ISari fc Esinl l200ll . A9), where ot is the Thomson scattering cross section, R is the 
shock radius, n is the external circumburst density in units of cm -3 , and r = axRn/3 is the 
Thomson optical depth of the radiation region 0. Another form of this ratio which has been 
used is 

F2 = fCjfZx = [4(p - 1)/ (p - 2)]t , (6) 



( IKobayashi et al.ll2007l ). It is apparent that Fl is smaller than F2 by a factor [4(p— 1) / (p— 2)], 
which can give substantial differences in analytical estimates of the IC spectral flux at its 
peak. It is worthwhile therefore to clarify the reason for the discrepancy. 

The Fl form (equation [5]) of the peak flux ratio is derived from an integral over the 
electron energy distributio n and a power-la w seed synchrotro n spectrum (See Eqn. (7.28), 
Rvbicki fc LightmaiJ[l979l : Also Eqn. (Al). ISari fc Esinl bpOlh . 



The F2 form (equation [B]) of the peak flux ratio, on the other hand, is obtained by 
solving for f^axl fmax f rom an equation ([7]) which relates the Compton Y parameter to the 
ratio of the luminositi es produced by the first order IC and the synchrotron mechanisms 
kobavashi et alil2007h . 



y - L IC ,ist/ L syn ~ v peak 



!C fv C {^ak)l V veakfv{y V eak) =2rlT lml , 



(7) 



3 The optical depth here differs from the usual definition by a factor 1/3, but for consistency with the 
usage in the literature (jPanaitescu fc Kumarll2000l ; IKobayashi et al.ll2007l) we keep the factor 1/3 here. 
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where vJ^ ak and v pea k are the peak frequency of the IC and synchrotron spectra, respectively, 
and fl c {^pe a k) and f v {y vea k) are the peak fluxes corresponding to the IC and synchrotron peak 
frequencies. The corresponding analytical approximation to the IC spectrum is a simpler 
one than in the previous Fl case, in that it is a broken power law (without the logarithmic 
corrections). In the low- frequency part of the IC spectrum, the broken power law is a good 
approximation to a numerically calculated IC spectrum. However, in the high-frequency part, 
the broken power-law analytical approximation under-estimates the numerical IC spectrum, 
which is larger (and has a flatter spectral index) than the broken power-law prediction (this 
under-estimation is avoided in the Fl case by including the logarithmic correction). In 
the F2 pure broken power law case, therefore, in order to keep the frequency integrated 
total luminosity radiated by the IC mechanism equal to the numerically computed total 
IC luminosity, and to preserve the desirable simple broken power law shape, an artificially 
boosted peak flux ratio is adopted, which leads to the same IC-to-synchrotron luminosity 
ratio. Thus, while the total energetics are the same for both analytical approximations, 
the IC flux expected over the GLAST (and AGILE) energy range differ. Whereas Fl is 
simpler for quick estimates since it involves pure power laws and correctly describes the 
global energetics, Fl with the logarithmic corrections to the power laws is preferable for 
more accurate GeV spectral flux estimates. 



